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ABSTRACT 

The theory of moving reversing layers for hot stars is updated to include an extensive line list, a radiative boundary 
condition from static model atmospheres, line transfer by scattering, and continuation to supersonic velocities. A Monte 
Carlo technique determines the theory's eigenvalue J, the mass flux, and the derived J's are in good agreement with 
the wind models of Pauldrach et al. (2001). The solutions' sensitivity to the photospheric microturbulent velocity Vt 
reveals that this parameter has a throttling effect on J: turbulent line-broadening in the quasi-static layers reduces 
the radiation force available to accelerate matter through the sonic point. If photospheric turbulence approaches sonic 
velocities, this mechanism reduces mass loss rates by factors ^ 3, which would partly account for the reduced rates 
found observationally for clumpy winds. 
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1. Introduction 

The line-driven winds of hot stars are unstable (Lucy & 
Solomon 1970, [LS70]; Owocki et al. 1988) and the result- 
ing density inhomogeneities strengthen emission lines such 
as Ha and Hen A4686, thus positively biasing empirical es- 
timates of the mass-loss rate (Lucy 1975; Abbott et al. 
1981). This bias could be corrected if the statistical prop- 
erties of the inhomogeneities could be predicted, but this is 
still not feasible. In consequence, dumpiness is now an es- 
sential extra parameter in diagnostic analyses of the super- 
sonic outflows from early-type stars (Schmutz 1997; Hillier 
& Miller 1999). That such investigations confirm dumpi- 
ness is no great surprise in view of work in the 1980's con- 
cerning X-ray emission, the black troughs of saturated P 
Cygni lines, and the variability of line profiles. But the 
magnitude of the concomitant reductions of the inferred 
$'s is a major surprise, as is the implication that theory 
overpredicts $'s by substantial factors. 

The "J's of Galactic O stars are investigated in the recent 
papers of Bouret et al. (2005) and Fullcrton et al. (2006), 
both of which also review earlier literature. Bouret et al. 
analyze the spectra of two 04 stars and derive <&'s that are 
factors of 3 and 7 below predicted values. On this basis, 
they conclude that our understanding of mass loss by O 
stars is in need of fundamental revision. 

Fullerton et al. study a much larger sample of 40 stars 
but, instead of analysing entire spectra, they focus on the 
Pv resonance doublet which, being seldom saturated, more 
readily yields reliable <E>'s than do lines from abundant 
species. Strong clumping, reducing previous empirical $'s 
by factors of 10 or more, is required to match the rates de- 
rived from P V. These authors therefore conclude that the 
"standard model" for hot-star winds needs to be amended. 

The theoretical <3?'s tested by Bouret et al. (2005) are 
those derived by Vink et al. (2000) with a refined version 
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of the Monte Carlo (MC) technique of Abbott & Lucy 
(1985, [AL85]). As such, these $'s are not eigenvalues of 
the differential equations governing outflows but are semi- 
empirical estimates based on an assumed wind structure. 
Accordingly, if these $'s are in conflict with observation, 
the fault may lie in the adopted stratification rather than 
in the underlying theory. 

At present, the closest to a complete theory of 
radiatively-driven winds is embodied in the code WM- basic 
of Pauldrach et al. (2001). Thus, for a critical test of wind 
theory, the above reduced empirical <E>'s should be com- 
pared to the eigenvalues $ of WM- basic models tailored 
to individual stars. Pending this, some progress is possi- 
ble with the models in Table 5 of Pauldrach et al. (2001). 
Interpolation using the Bouret et al. (2005) parameters for 
HD 190429A gives a WM-basic $ ~ 6.3 x 1O~ 6 yr' 1 , a fac- 
tor 3.5 larger than their clumped- wind model. Thus claims 
that line-driven wind theory is in need of major revision 
appear to be well founded. 

Recent papers on dumpiness confirm what was already 
evident in the mid 1980's, namely that the "standard 
model" of line-driven winds must eventually be replaced by 
3-D time-dependent gas-dynamical models with simultane- 
ous NLTE radiative transfer. Such models are necessary if 
we wish to predict from first principles spectral features 
formed at supersonic velocities. But if clumping is negligi- 
ble in the neighbourhood of the sonic point, the standard 
model's predictions of $ should be fairly accurate since the 
differential effect of dumpiness at supersonic velocities on 
the radiation field at the sonic point is surely slight. Thus 
the evidence that the predicted <&'s are far from accurate 
perhaps indicates that dumpiness is already significant at 
the sonic point, a conclusion with some observational sup- 
port (Bouret et al. 2005). 

Notwithstanding this latter possibility, this paper 
is motivated by the belief that the eigenvalues 4> of 
stationary-wind solutions should indeed be fairly accurate. 
Accordingly, the moving reversing layer model of LS70 is 
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revived in order to compute mass fluxes for star atmo- 
spheres and to investigate their sensitivity to assumptions 
and parameters. 



2. An outflowing reversing layer 

In this section, the LS70 model of stationary flow acceler- 
ating to sonic velocity is updated. The major changes are a 
vast increase in the line list, a realistic lower boundary con- 
dition on the radiation field, an improved treatment of line 
formation, and continuation to supersonic velocities. The 
solution procedure is also changed. Because line formation 
by pure absorption is replaced by scattering, solution with 
initial-value integrations is no longer possible. Accordingly, 
a MC technique akin to that of AL85 is developed. 

Despite these improvements, the resulting code falls far 
short of the physical detail and numerical precision achieved 
by codes such as WM- basic, TLUSTY (Lanz & Hubeny 
2003) , CMFGEN (Hillier & Miller 1998), and FASTWIND 
(Puis et al. 2005). But high accuracy is surely not needed in 
identifying an effect capable of reducing $ by factors <; 3. 



2.1. The model 

As in LS70, a Schuster-Schwarzschild model is adopted. 
Thus the lower boundary is a continuum-emitting surface 
above which is an outflowing plane-parallel isothermal layer 
where only line- and electron scatterings occur. In addition 
to chemical composition, the basic parameters are effective 
temperature T e ff, surface gravity g, microturbulcnt veloc- 
ity vt, and mass flux J = $/47ri? 2 . 

The composition is solar with Nhc/Nh = 0.1 (Grevess 
& Sauval 1998), with the included ions as in Table 1 of 
Lanz & Hubeny (2003). The temperature of the isothermal 
flow is set = 0.75T e ff, and H and He are assumed to be 
fully ionized, so that thermodynamic coefficients are height- 
independent. 



2.2. Gas dynamics 

The structure of the reversing layer is determined by 
the equations governing stationary isothermal flow. With 
sphericity neglected, the continuity equation has the inte- 
gral 

pv = J (1) 



and the equation of motion can then be written as 
2 2 1 dv 



(2) 



Here a = y/P/p is the isothermal speed of sound and the 
effective gravity 



9eff =g-g e 



gi 



(3) 



where g e and gi are the radiative accelerations due to 
electron- and line scatterings, respectively. 

From Eq.(2), we see that the sonic point v = a is a 
critical point at which the velocity gradient is finite only 
if g e ff = 0. This condition was used in LS70 to determine 
the eigenvalue J. 

Subsequent to LS70, Castor, Abbott & Klein (1975, 
[CAK]) pioneered the development of a complete theory 



of stationary line-driven winds by exploiting the compu- 
tational economy of the Sobolev approximation. However, 
the critical point is then no longer at the sonic point but 
is displaced downstream into the supersonic flow. Here, as 
in LS70, the Sobolev approximation is not employed. In 
consequence, gi is a functional of the solution, and so the 
critical point remains at the sonic point (Lucy 1975; Poe et 
al. 1990). 

2.3. Stratification 

In AL85 and Vink et al. (2000), the MC transfer is per- 
formed in an observationally-motivated kinematic model, 
with $ determined iteratively by imposing a global dynam- 
ical constraint. Here, with the aim of achieving greater lo- 
cal consistency, a family of dynamical models is generated 
corresponding to an assumed functional form for ge(v). The 
mass flux J is then determined iteratively by bringing ge(v) 
and its MC estimator gi{v) into approximate agreement. 
The adopted ge(v) is the two-parameter formula 



= g*max 



with g„ = g - g e 



(4) 



With < S < 1 and s > 0, this formula is such that 
g e ff — at v = a and is negative at supersonic velocities. 
The discontinuous derivative of this function at the join 
of its two segments allows the sharp increase in gt as the 
outflow accelerates up to and through the sonic point to be 
captured. 

The solution of Eq. (2) that avoids a singularity at the 
sonic point has 



dv 

! 

dx 



1 / dlng\ 1 



Accordingly, the velocity law v(x; 5, s) is derived from Eq. 
(2) by inward and outward integrations from v = a, each 
with initial derivative given by Eq. (5). 

Together with Eq. (1) and the isothermal assumption, 
v(x; 5, s) determines the stratification throughout the com- 
putational domain (xo,Xi). The lower boundary x is at 
electron-scattering optical depth r e = 2/3 below the sonic 
point. The upper boundary x\ is at v ~ 4a. 

Note that, as in solar-wind theory, the solution has an 
X-type topology at the critical point. This is assumed with- 
out proof because of the difficulty of analysing the solution 
topology when gi includes contributions from scattered ra- 
diation (Poe et al. 1990). 

2.4. Radiation field 

The frequency distribution of the continuum intensity J+ 
emitted at xq is that of the emergent continuum flux of the 
TLUSTY (Lanz & Hubeny 2003) solar-abundance model 
of the specified T e ff and g; the ionization edges due to H 
and He are therefore accurately represented. The angular 
distribution is assumed uniform. 

At xi, the integrated flux of the outwardly-directed ra- 
diation I„{p) is constrained to = crT^ff. Correspondingly, 
I~ = at xi, so that back-scattering from the highly su- 
personic flow is neglected. Nevertheless, x\ is far enough 
downstream that, for v ~ a, each line's self-interactions 
with back-scatterings from v > a are treated. 
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As it propagates through the reversing layer, the contin- 
uum is modified by electron scatterings and line transitions. 
Electron scatterings are assumed to be isotropic but fre- 
quency shifts in the co-moving frame (cmf) due to the elec- 
trons' thermal motions are taken into account. Following 
Sect. Ilf in AL85, line transitions are treated as resonance 
scatterings. These are assumed to be isotropic with a re- 
emitted frequency in the cmf in accordance with complete 
redistribution. The line profile <f> v comprises a Doppler core 
due to thermal and microturbulent motions and wings due 
to radiation damping. 

2.5. Ionization and excitation 

The treatment of ionization and excitation is designed to 
give moderately accurate populations for the lower levels 
of transitions that are major contributors to ge- Thus the 
emphasis is on the ground states and low-lying metastable 
levels of abundant ions. Conveniently, NLTE departure co- 
efficients for O-star atmospheres can be downloaded from 
the TLUSTY website, and data from a model with T e f f — 
50000if , logg — 4 are plotted in Fig. 1 for the important 
ion Fe v. This diagram shows the expected convergence to 
rigorous LTE in deep layers and the growing importance of 
NLTE in the surface layers. However, because of collisional 
coupling, the low-lying superlevels of Fe v remain approxi- 
mately in LTE relative to the ground state. In contrast, for 
the ionization ratio Fe v/Fe iv, significant departures from 
LTE occur for n e ^ 13 dcx. Nevertheless, changes in n e 
remain the dominant cause of ionization gradients in the 
surface layers. 

In view of the above, excited levels are assumed to be in 
LTE with their ground states, whose departure coefficients 
are fixed at the TLUSTY values in a representative layer. 
Thus, the ionization equation is 

n J+1 n e _ / &i„/ + i \ / n J+1 n e \ * 
n J \ h. j / t V n J J f 

where * denotes LTE and f indicates values for the layer 
with T e ~ 0.75T e// . 

2.6. Line list 

The basic line data is from Kurucz & Bell (1995). A typ- 
ical model treats the line formation of ~ 10 5 lines in the 
wavelength interval 228 - 4000A . 



3. Monte Carlo code 

In this section, a MC code is described for calculating ge as 
a function of height in the outflowing reversing layer defined 
in Sect. 2. The description is brief, focussing on aspects not 
covered in AL85. 

3.1. Discretization 

The continuous model derived in Sect. 2 is approximated 
by ~ 100 slabs, within each of which p and v are constant. 
The velocity increment between neighbouring slabs is ev- 
erywhere <C Vt- 
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Fig. 1. Ratios of departure coefficients relative to the ground 
state for selected levels of Fe as functions of n e for TLUSTY 
model G50000g400vl0.12. Ratios are plotted for superlevels of 
Fe v with excitation energies < 10 eV. Also plotted is the ratio 
of b- values for the ground states of Fe v and Fe IV. The filled 
circle is the point where T e ~ 0.75T e //. 

3.2. Scattering histories 

A total of N photon packets of energy eo are launched into 
the computational domain at xq in time At. Their initial 
rest frame frequencies v sample the continuum flux distri- 
bution of the TLUSTY model and their initial direction 
cosines are given by p — y/z, where here and below z de- 
notes a random uniform variate obtained with ran2 (Press 
et al. 1992). 

A packet's straight-line trajectory in a slab terminates 
after travelling distance s = min(sb, s e , si), where Sb is 
the distance to the slab's edge, s e = —£n(z)/a e , and 
se = —£n(z)/£ v i. Here <7 e and l v * are, respectively, the 
electron- and line scattering coefficients per unit volume, 
the latter evaluated at v' , the packet's cmf frequency. Note 
that these coefficients are constant within each slab. 

If the packet undergoes an electron scattering, its new 
p = 1 — 2z and its new v' is obtained by sampling the angle- 
averaged redistribution function (Mihalas 1978, p. 432). For 
a line scattering, the new p = 1 — 2z and the new v' is 
obtained from the assumption of complete redistribution. 

The scattering history of a packet is complete when it 
exits the computational domain at xq or x\. The number 
J\f + escaping at x\ constitute the emergent radiation and so 
must have an integrated flux corresponding to T e ff. Thus, 
at the completion of the MC experiment, the packets' en- 
ergy e is given by 

^^=^ft (7) 

3.3. Estimator for ge 

With e a /At from Eq. (7), the MC radiation field, compris- 
ing the trajectories of Af packets, can be transformed to 
physical units and then used to calculate ge. 
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If the intensity I v (x,n) were known, gg could be com- 
puted from the formula 



P9t = — 
c 



i v 'I v pL d\i dv 



(8) 



1 Jo 



valid for v <C c. This can be applied to the energy-packet 
model using the intensity estimator 



I v dv dbJ = 4— T7 ^ — s 

At V f-f e 



(9) 



where the summation is over all pathlength segments s in 
the reference volume V of packets in (y, v + dv) propagating 
in solid angle element dtu. This formula is derived in Lucy 
(2005) and repeated in order to correct a typographical 
error. Eq. (9) relates the energy-packet and the specific- 
intensity representations of radiation. 

In the present application, v <C c and so setting e„ = eo 
is an excellent approximation. Also the geometry is plane- 
parallel, so V is a column of unit cross-section and length 
Ax, the slab's width. Substitution into Eq. (8) then gives 
the estimator 



P9t = 



1 e 1 

c At Ax 



tv' P ■ 



(10) 



where the summation is over all trajectories in the slab. 

This optimum estimator for gg exploits all the infor- 
mation in the packets' trajectories. This would not be the 
case if gg were derived from the momentum transfers oc- 
curring when packets undergo line scatterings. Even if no 
line scatterings occur in a slab, Eq.(10) returns a non-zero 
gg, provided only that at least one packet propagates with 
a v' at which l v i ^ 0. 

Because this calculation of gg does not rely on the 
Sobolev approximation, the effects of scattered radiation 
discussed by Owocki & Puis (1999) are included. Moreover, 
line-profile overlapping and multiline transfer are accu- 
rately treated. These effects are excluded when Sobolev 
theory is used in the narrow-line limit with the assump- 
tion that each line interacts with unattenuated continuum. 



3.4. Computational economy 

With the procedure as described, a prohibitive amount 
of computer time would be required for an accurate gg. 
Accordingly, two simplifications are introduced. 

A MC treatment of re-emission with complete redistri- 
bution should first randomly select which of the overlapping 
lines contributing to t v i is the actual absorber and then 
sample its profile <j> v i to select the v' of the emitted packet. 
This has been simplified by always taking the absorber to 
be the largest contributor to i v i, which is identified when 
l v i is tabulated prior to the MC calculation. Little loss of 
accuracy is expected since absorption usually occurs in a 
line's Doppler core and this typically strongly dominates 
the damping wings of neighbouring transitions. 

The second simplification avoids trapping in the 
Doppler cores of strong lines. A photon undergoing numer- 
ous scatterings in a Doppler core is described by Rybicki & 
Hummer (1969) as being imprisoned at a point and from 
which it escapes only when the low probability event of 
emission in the damping wings occurs. When a MC packet 
is similarly imprisoned, its escape via the damping wings 



is expedited. The actual procedure is as follows: with stan- 
dard notation, the Voigt profile is approximated (Mihalas 
1978, p. 281) by the gaussian core when \v\ < v* and by the 
damping wings when |u| > v*, where v* is such that the 
two terms are equal. If a slab is optically thick in a Doppler 
core, the next re-emission is forced to occur in the wings. 
Since <p v cx 1/v 2 in the wings, a random v = ±v*/z. 



4. Examples 

In this section, moving reversing layers are derived for 
the three TLUSTY atmospheres with parameters match- 
ing those of models in Table 5 of Pauldrach et al. (2001). 
Mass fluxes can then be compared without interpolation, 
thereby testing the approximations used to construct the 
simple code described in Sects. 2 and 3. 



4.1. Parameter fitting 

With T e f f , g and composition fixed, the remaining parame- 
ters are J, 5, s and v t . In the WM-basic models, v t increases 
from 10 km s^ 1 at the photosphere to 0.1 in the terminal 
flow. Here, we ignore this gradient and set Vt = 10 km s _1 . 

With a trial parameter vector (J, 6, s), a MC calculation 
is carried out as described in Sect. 3. Eq.(10) then gives gg 
throughout the reversing layer and thus allows us to see how 
closely the assumed gg(v) given by Eq. (4) is recovered. The 
goodncss-of-fit is quantified by calculating 



9i ~ 9e 



9i 



(11) 



where the summation is over the / slabs with 0.25 < v/a < 
1.75. Velocities outside this interval are excluded to avoid 
the edge effects due to imperfect boundary conditions. 

With many repetitions of this procedure, the vector 
(J, 5, s) minimizing r\ can be found. In the initial search, 
models are calculated with M = 10 , but this is increased to 
4x 10 7 in the close neighbourhood of the minimum and then 
further increased to M ^ 10 8 for the illustrations herein. 

Since the interval (0.25, 1.75) includes the sonic point, 
the regularity condition at v — a is satisfied exactly if rj = 0. 
But because of MC sampling errors and the limited flexibil- 
ity of the formula for gg, this is not achievable. Accordingly, 
minimizing r\ replaces the condition g e ff = at v = a as 
the criterion determining the eigenvalue J. 

4.2. Model D-50 

Model D-50 of Pauldrach et al. (2001) has T eff = 50000K, 
log 5 = 4, R/Rq = 12, and $ = 5.6 x 10~ 6 M Q /yr, corre- 
sponding to a photospheric mass flux J(gm s _1 ) = —4.40 
dex . The matching TLUSTY model is G50000g400vl0.12, 
details of which were used for Fig. 1. From the emergent 
spectrum of this model, continuum fluxes were estimated 
at intervals of 0.1 in logv and at both sides of ionization 
edges. Linear logarithmic interpolation then defines the v- 
distribution of at xq. 

The best-fit moving reversing layer has parameter vec- 
tor (log J, 6, s) = (-4.41,0.61,2.03) which gives r\ = 
0.048, and this fit worsens with displacement vectors 
(±0.05, ±0.03, ±0.05). Thus the mass flux of this model 
agrees with the WM-basic model to within 0.05 dex. 
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Fig. 2. Comparison of radiative accelerations due to lines. The 
filled circles are values given by the MC estimator gt, and the 
solid line is the function gi(v) assumed in deriving the reversing 
layer's stratification. The open circle indicates the sonic point. 



The success in reproducing the assumed gt is shown in 
Fig. 2. The agreement is seen to be excellent from v ~ 10 to 
~ 70 km s _1 , an interval containing the sonic point at 22.6 
km s . In particular, the discontinuity in the derivative of 
gt(v) at v = 17.5 km s _1 closely matches a sharp change 
in the slope of gt(v). On the other hand, the fit is less 
successful for v & 10 km s _1 , where Eq. (4) reduces to 
gt =g*6. 

4.3. Model D-40 

Model D-40 of Pauldrach et al. has T eff = 40000K, logg = 
3.75, R/Rq = 10, and $ = 0.24 x lO~ 6 M /yr, correspond- 
ing to J(gm s _1 ) = —5.60 dex . In this case, the source of 
b-values at the representative point and of the continuum- 
flux distribution is the TLUSTY model G40000g375vl0.12. 

The best- fit moving reversing layer for v t — 10 km s _1 
has (log J, 6, s) = (-5.65,0.58,2.55), with r\ = 0.048. Thus 
the mass flux is lower than that of the WM- basic model by 
a barely significant 0.05 dex. 

4.4. Model S-30 

Model S-30 has T eff = 30000K, log g = 3.00, R/Rq = 
27, and <f> = 5.0 x 10~ 6 Mq/ut, corresponding to 
Jigms- 1 ) = -5.15 dex . The matching TLUSTY model is 
G30000g300vl0.12. 

The best- fit moving reversing layer for vt = 10 km s _1 
has (log J, 5, s) = (-5.45,0.79,1.75), with r] = 0.072. Thus 
the mass flux is lower than that of the WM-basic model by 
an apparently significant 0.3 dex. However, the reliability 
of this model has been questioned (Puis et al. 2005, p686). 



4.5. Consistency checks 

The good agreement of the J's with those of the reliable 
WM-basic models provides strong support for the simpli- 
fying assumptions introduced in Sects. 2 & 3. A further 
consistency check is obtained from the lines contributing 
to gt. The treatment of ionization and excitation in Sect. 
2.5 focuses on low-lying energy levels in the expectation 
that the bulk of the radiative driving derives from scatter- 
ing by such levels. This is tested for the layer with v ~ a 
in model D-50. For this layer, gt ~ g*, with the fractions 
0.70 and 0.88 arising from levels with excitation energies 
< 5 and < 10 eV, respectively. 

5. Variation of parameters 

In this section, insight is sought into the behaviour of mov- 
ing reversing layers by varying parameters. 

5.1. Mass flux 

The excellent fit of gt to gt in Fig. 2 is lost if parameters 
change. In Fig. 3, the effects of the displacements AJ = 
±0.3 dex are shown. The changes in gt are in the expected 
sense. Thus an increase in J results in a gt that falls below 
the assumed gt(v) and is therefore incapable of accelerating 
the outflow. 

Fig. 3 supports the assumption of an X-type topology 
since there is no hint of the solution indeterminacy char- 
acteristic of a nodal-type topology (Poe et al 1990). The 
sensitivity to mass flux is consistent with there being a 
unique eigenvalue J for which there is a stationary outflow 
accelerating from sub- to supersonic velocities. 

5.2. Microturbulence 

In a second sensitivity experiment, the effect of varying vt 
is investigated. Fig. 4 plots gt(v) for v t — 6.7 and 15 km 
s _1 . In this case, the results are less readily anticipated. 

For w < 15 km s _1 , gt increases with w t . This is a con- 
sequence of the broadened lines in the quasi-static layers 
intercepting an increased fraction of the continuum. A sim- 
ple analysis of this effect is given in Sect. lib of LS70. With 
a slowly-varying logarithmic factor neglected, this analy- 
sis shows that, in a static atmosphere with a fixed col- 
umn density of absorbing ions, a strong resonance line has 
gt oc Afjj, which is cx v t when microturbulence dominates 
the ions' thermal motions. 

Following a crossover at v ~ 20 km s _1 , this depen- 
dency is reversed and gt becomes an decreasing function of 
Vt- This is again a consequence of the broadening of lines 
formed in the quasi-static layers. Continuum radiation that 
would have been available to accelerate matter moving with 
v ~ 30 — 40 km s _1 has already been partly absorbed in 
the turbulent flow at subsonic bulk velocities. 

The final feature to be explained is the convergence of 
the gt plots for v ^ 60 km s _1 . This is the Sobolev regime. 
Matter has been accelerated from out of the shadow of lines 
formed in the quasi-static flow and the narrow-line limit is 
approached in which gt becomes independent of (j) u and 
therefore of v t . 

With Fig. 4 understood, the next point of interest is 
the change in (J, d, s) required to recover a good fit be- 
tween gi and gt. The above discussion indicates that as v% 
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Table 1. Models with varying microturbulence 



Velocity (km/s) 

Fig. 3. Effect of varying mass flux J, with other parameters 
as in Fig. 2. The filled and open circles are values of ge when 
J(gm s _1 ) = —4.11 and —4.71 dex, respectively The solid line 
is the function ge(v) assumed in deriving the reversing layer's 
stratification, and the large open circle indicates the sonic point. 



increases so will S, as lines then provide more support for 
the quasi-static layers, but J will decrease to compensate 
for the reduced flux available to accelerate matter through 
the sonic point. 

For each of the three WM-basic models considered in 
Sect. 4, a sequence of best-fit models with varying v t is 
given in Table 1. Each sequence has the continuum flux 
distribution and the b-values of the TLUSTY model iden- 
tified in Sects. 4.2-4.4. In consequence, there is an inconsis- 
tency since the TLUSTY models assume v t = 10 km s _1 . 
But the dominant effects of v t on the parameter vector are 
surely those just discussed. Note that the effect of v t on 
line-blocking and hence backwarming is approximately ac- 
counted for by Eq. (7) . Thus, because packets that re-enter 
the photosphere are implicitly thermalized and re-radiated, 
the continuum flux at x exceeds ^T^ff by an amount that 
depends on v t . 

From Table 1, we see that S and J vary in the sense 
expected. In particular, J decreases by ~ 0.3dex as v t in- 
creases by 0.18dex, corresponding to a power-law exponent 
of —1.7. However, although this effect is readily understood 
and undoubtedly real, the quality of the fits is markedly 
poorer for the D-50 models with v t = 15 and 22.5 km s _1 . 
This indicates the need for a more complicated fitting func- 
tion gi(v) to replace Eq. (4). 

Turbulent broadening is treated here in the microtur- 
bulcnt limit in which the length scale of random velocity 
fluctuations is -C the mean free paths of line photons. But 
the mechanism does not require this assumption: any pho- 
tospheric line-broadening that reduces pgi at v ~ a will 
contribute to this throttling effect. In particular, therefore, 
turbulence with length scales up to ~ the photospheric scale 
height might well be important, and there is strong obser- 
vational evidence (Conti & Ebbets 1977; Howarth et al. 
1997) for turbulent velocities well in excess of the standard 
Vt — 10 km s^ 1 required by line strengths. 
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Fig. 4. Effect of varying microturbulent velocity vt, with other 
parameters as in Fig. 2. The filled and open circles are values 
of gt when v t = 15 and 6.7 km s _1 , respectively The solid line 
is the function ge{v) assumed in deriving the reversing layer's 
stratification, and the large open circle indicates the sonic point. 



6. Discussion 

The above calculations reveal a fundamental coupling be- 
tween line formation in the quasi-static photospheric layers 
and the mass flux J that can be accelerated to supersonic 
velocities. This low- velocity, non-Sobolev effect must be ac- 
curately treated when J or $ arc computed as eigenvalues 
of the equations of stationary line-driven winds. 



6.1. Previous work 

This coupling was already recognized in LS70. In Fig. 2 of 
that paper, the precipitous drop of J at logT e ff ~ 4.5 was 
attributed to the absorption of radiation at i^(l + a/c) by 
the damping wings of the Civ doublet. Here, with the inclu- 
sion of numerous weak lines, J falls less precipitously when 
photospheric line broadening reduces the contributions of 
strong lines to pgg at v ~ a. Moreover, rather than damp- 
ing wings, turbulent broadening of Doppler cores is a more 
plausible cause of a reduced radiation force in the transonic 
flow. 

Because the throttling effect of turbulent broadening 
requires non-Sobolev line transfer, it is absent from CAK 
theory. The ion thermal speed v t h that appears in CAK 
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formulae is for dimensional reasons and scales out when 
derived properties are expressed in physical units (Abbott 
1982, Poe et al. 1990). 

The hydrodynamic structure of WM- basic models are 
obtained with CAK theory and thus a dependence of the 
eigenvalue <!> on photosphcric turbulence would seem to be 
excluded. However, the CAK force multiplier parameters 
are iteratively adjusted to fit gi given by a detailed NLTE 
transfer code that assumes a velocity-dependent microtur- 
bulcnce with w t = 10 km s _1 at the photosphere. Although 
this introduces a dependence of $ on the displacement 
of the CAK critical point into the supersonic flow suggests 
that the throttling effect will be less powerful than found 
here. But if the CAK formalism were eliminated and the 
stratification derived directly from the gas dynamic equa- 
tions with critical point at the sonic point, an accurate 
treatment of the dependence of $ on v t would be possi- 
ble with this code. 

Although not identified with turbulence, the throttling 
effect of an enhanced v t h/a is evident in the pure absorp- 
tion calculations reported by Poe et al. (1990). Moreover, 
they summarize the implications of this parameter for time- 
dependent simulations and the uniqueness of stationary so- 
lutions. 

Hillier et al. (2003) identify the treatment of microtur- 
bulence in calculating pgi near the sonic point as a major 
unsolved issue and recognize that the Poe et al. work can be 
re-interpreted to anticipate the effects of microturbulence. 

6.2. Calculating $ 

Two theoretical techniques for computing $ have emerged, 
one fundamental, the other semi-empirical. The fundamen- 
tal approach determines $ or J as an eigenvalue of the 
governing differential equations by demanding regularity 
at the critical point. The semi-empirical approach, on the 
other hand, determines <f> algebraically by demanding that 
the flux of mechanical energy in the terminal flow should 
be accounted for by the loss of radiative luminosity due to 
line- and electron scatterings in the outflowing matter. 

A consequence of this investigation is that the funda- 
mental approach is harder than hitherto realized. An es- 
sential phenomenon, photospheric turbulence, is not pre- 
dictable as yet from first principles and only measurable 
with moderate accuracy, even in the favourable but unlikely 
circumstances that such motions are restricted to the micro- 
turbulent limit and do not vary with height. Accordingly, 
by default, the semi-empirical technique is currently prefer- 
able as a source of $'s for input into other astrophysical 
investigations. The details of transonic flow are then incon- 
sequential because the kinetic energy density at the sonic 
point is negligible compared to that of the terminal flow 
(w^/2) which, moreover, is an observable. 

As discussed in Sect. 1, the semi-empirical $'s obtained 
by Vink et al. (2000) are too high according to Bouret et 
al. (2005). This may well be because the AL85 assump- 
tions of laminar flow with a monotonic velocity law are re- 
tained. But these assumptions are not fundamental to the 
MC technique. The Vink et al. (2000) investigation could 
usefully be repeated with the stratification updated in ac- 
cordance with the already cited papers on dumpiness as 
well as the recent papers of Puis et al. (2006) and Oskinova 
et al. (2006). Moreover, the treatment of line formation 
could be improved to treat branching (Sim 2004). 



A survey of $'s for hot stars could also be carried out 
with the diagnostic code CMFGEN, which already incor- 
porates dumpiness and random non-monoticity. For radii 
where the bulk velocity v(r) is highly supersonic, we have 
to a good approximation 

\v 2 = J r {9i - 9*)dr (12) 

By demanding that this equation is satisfied at r = oo, $ 
could be determined algebraically in a manner analogous to 
AL85 and Vink et al. (2000). Further wind parameters, such 
as (3 in the velocity law, could be determined by minimizing 
violations of this equation over its range of validity. 

7. Conclusion 

The aim of this paper has been to examine claims that the 
theory of line-driven winds overpredicts mass loss rates. To 
this end, the theory of moving reversing layers has been 
updated and used to investigate the sensitivity of transonic 
flow to input parameters. The reported calculations reveal 
that photosphcric turbulence, by reducing the driving force 
at the sonic point, plays a decisive role in regulating the 
mass flux of stationary flow. Codes that compute $ as 
an eigenvalue will therefore overpredict if turbulent line- 
broadening is underestimated. This may occur as a result 
of observational underestimates of vt or because the mi- 
croturbulent limit is violated. A further possible cause of 
overprediction is reliance on the CAK critical point in deter- 
mining the eigenvalue. Because this point is at supersonic 
velocity, the throttling effect of photospheric turbulence is 
likely to be reduced, a topic meriting further research. 

The calculations presented in Table 1 suggest the con- 
jecture that the mass loss rates of hot stars are less than 
predicted because of the throttling effect of photosphcric 
turbulence. With the critical point at the sonic point as re- 
quired by non-Sobolev line transfer, this effect reduces J for 
stationary flow by a factor <; 3 for near sonic turbulence. 

However, additional work is required to establish this 
conjecture. The MC technique presented herein is barely 
adequate for computing J's when vt 15 km s . 
Repetition with an improved parameterization of gi or so- 
lution with an alternative procedure is therefore desirable. 
But, more fundamentally, the conjecture is motivated by 
the behaviour of stationary solutions, and these are unsta- 
ble (Sect.l). Ultimately, therefore, as the referee remarks, 
time-dependent calculations are needed to see if the time- 
averaged J is similarly constrained by photosperic turbu- 
lence. 

With regard to analyses of individual stars, the eigen- 
value $ = 4irR 2 J obtained from a moving reversing layer 
computed with the spectroscopic v t should be compared 
with the empirical $ obtained from the supersonic wind 
with dumpiness and non-monoticity included. Such anal- 
yses could usefully exploit the wide range of Vt's reported 
for O stars (Bouret et al. 2003; Heap et al. 2006). 

Acknowledgements. I am grateful to the referee, J. Puis, for detailed 
comments on the original version. 

References 

Abbott, D.C. 1982,ApJ, 259, 282 



8 



Lucy: Mass fluxes 



Abbott, D.C., Bieging, J.H., & Churchwell, E. 1981, ApJ, 250, 645 

Abbott, D.C., & Lucy, L.B. 1985, ApJ, 288, 679 (AL85) 

Bouret, J.-C, Lanz, T., Hillier, D. J., Heap, S. R., Hubeny, I., Lennon, 

D. J., Smith, L. J., & Evans, C. J. 2003, ApJ, 595, 1182 
Bouret, J.-C, Lanz, T., & Hillier, D. J. 2005, A&A, 438, 301 
Castor, J.I., Abbott, D.C., & Klein, R.I. 1975, ApJ, 195, 157 (CAK) 
Conti, P. S., & Ebbets, D. 1977, ApJ, 213, 438 
Fullcrton, A.W., Massa, D.L., & Prinja, R.K. 2006, ApJ, 637, 1025 
Grevess, N., & Sauval, A.J. 1998, Sp. Sci. Rev., 85,161 
Heap, S. R., Lanz, T., & Hubeny, I. 2006, ApJ, 638, 409 
Hillier, D.J., & Miller, D.L. 1998, ApJ, 496, 407 
Hillier, D.J., & Miller, D.L. 1999, ApJ, 519, 345 

Hillier, D.J., Lanz, T., Heap, S. R., Hubeny, I., Smith, L. J., Evans, 
C. J., Lennon, D. J., & Bouret, J. C. 2003, ApJ,588,1039 

Howarth, I. D., Sicbcrt, K. W., Hussain, G. A. J., & Prinja, R. K. 
1997, MNRAS, 284, 265 

Rybicki, G. B., & Hummer, D. G. 1969, MNRAS, 144, 313 

Kurucz, R.L., & Bell, B. 1995, Kurucz CD-ROM No.23 

Lanz, T. & Hubeny, I. 2003, ApJS, 146, 417 

Lucy, L.B. & Solomon, P.M. 1970, ApJ, 159, 879 (LS70) 

Lucy, L.B. 1975, Mem.Soc.R.Sci.Liege,8,359 

Lucy, L.B. 2005, A&A, 429, 19 

Oskinova, L. M., Fcldmeier, A., & Hamann, W.-R. 2006, MNRAS, 
372, 313 

Owocki, S.P., Castor, J.I., & Rybicki, G.B. 1988, ApJ, 335, 914 
Owocki, S. P. & Puis, J. 1999, ApJ, 510, 355 

Pauldrach,A.W.A., Hoffmann,T.L., & Lennon,M. 2001, A&A, 375, 
161 

Poe,C.H., Owocki, S. P., & Castor.J.I. 1990, ApJ, 358, 199 

Press W.H., Teukolsky S.A., Vcttcrling W.T., Flannery B.P. 1992, 

Numerical Recipes. Cambridge Univ. Press, Cambridge 
Puis, J., Urbancja,M.A., Venero,R., Repolust,T., Springmann,U., 

Jokuthy,A., & Mokiem.M.R. 2005, A&A, 435, 669 
Puis, J., Markova,N., Scuderi,S., Stanghcllini,C, Taranova,O.C, 

Burnley,A.W., & Howarth,I.D. 2006, A&A, 454, 625 
Schmutz,W. 1997, A&A, 321, 268 
Sim, S. 2004, MNRAS, 349, 899 

Vink, J.S., de Koter, A., & Lamers, H.J.G.L.M. 2000, A&A, 362, 295 



